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We discuss the nature of pressure induced phase transitions in standard Quantum Paraelectrics 
near quantum critical point. From a microscopic theory we first show that near the critical point 
the transition temperature T c (p) varies as (f —p/p c )^ with f = |. Within the Landau scheme, it is 
shown that for £ < 1, a condition satisfied by most quantum paraelectrics the transition is inevitably 
first order in nature. The present analysis is valid for non-polar impurity induced transitions in these 
materials and is supported by some recent experiment. 
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I. INTRODUCTION AND SUMMARY 



II. MEAN FIELD ANALYSIS 
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The static dielectric susceptibility in the materials like 
SrTi03 and KTa03 docs not show any divergence down 
to zero Kelvin, it rather saturates at a very high value 
(O(10 4 )). Ferroelectric transitions in materials which 
are similar in structure, for example in BaTi03, are of 
displacive type and are well described by softening of a 
zone center transverse optical mode. This scenario, that 
is, absence of transition and high dielectric susceptibil- 
ity, has been attributed to smallness of the gap in the 
corresponding optical branch, smallness of the optical 
gap makes quantum fluctuations relevant in these ma- 
terials. These materials are, therefore, called quantum 
paraelectrics^. There is a revival of interest, particularly 
regarding the nature of phase transition at quantum crit- 
ical point in these materials^^. The effects of pressures^ 
and impurities on the dielectric susceptibilities have been 
well studied experimentally. One recent experimental^ 
report indicates that these materials show phase separa- 
tion near QCP. A general consensus arising form these 
experiments is that application of hydrostatic pressure 
moves them away from criticality and the possibility of 
transition is suppressed. One needs to apply some kind 
of negative pressure to induce phase transitions in these 
materials. One way to apply negative pressure is to 
put non-polar impurities which creates local pressure de- 
ficiencies and thus induce phase transitions. Remark- 
ably when experimentally^ one finds T c ~ (n — n c )2 
(where n is the average impurity concentration and n c 
is the critical value), the theoretically estimated^ tran- 
sition temperature for pressure induced transition goes 
as T c ~ (p + p c ) 5 (where p is hydrostatic pressure and 
p c is the critical value). This motivates us to develop 
a description, suitable for properties of pressure induced 
phase transition, which can be used to understand the 
newly found phase separation in ferroelectric transition 
near ferroelectric quantum critical point. In the follow- 
ing, using a microscopic theory in the paraclcctric phase, 
we first show that T c ~ (p + p c )^ within some reasonable 
approximation. Then we use this scaling behavior of T c 
in a Landau functional to discuss ferroelectric transition 
in such systems. 



Application of hydrostatic pressure (also the effect 
of non-polar impurity) will couple to optical mode via 
strain. In this case our starting Hamiltonian takes the 
form 



H 
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dq- 
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+ g dk dq e(k) u q u k - q - pi dqS(q)e(q) 
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Here terms in the parentheses describe optical phonons in 
harmonic approximation, the quartic term describes the 
anharmonic coupling between these phonons. Last three 
terms describe the result of application of pressure to 
lowest possible order. The strain field e(q) couples to unit 
cell displacement related to optic mode with amplitude 
<?, the harmonic acoustic phonons have force constant 
K and in the last term the coupling of the hydrostatic 
pressure p to the static strain has unit strength. The 
details of this hamiltonian are described in our earlier 
work^. On Integrating out the strain field, we get an 
effective Hamiltonian of the form 



H = dq 



\pI + \{ UJ t-v + gp + v5q 2 a 2 ) u q U- q 



-Xb J HidqiU qi u q2 Uq 3 U- qi - 



+ 4 A 



where the normalized quartic coupling 



A« = (A-f). 



(2) 
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Now with p q = ii q = —iLJUq in the kinetic energy term, 
and within the quasi harmonic approximation, for the 
quartic term, i.e. 
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the Hamiltonian is given by, 
H = \ ^{u 2 q -u 2 )u q u„ 



(5) 



Here uj q is renormalized optical frequency in the paraelec- 
tric phase where < u > is identically zero and is given 
by 



and 



u> 2 = u>q+ gp + vScosqa + 3\r<j 
~ lJ 1 — v + gp + vSa 2 q 2 + 3\r<t. 

o = ^(T U(? (0)u_,(0 + )) 



(6) 



(7) 



is the mean square fluctuations of displacement in optical 
mode. The susceptibility, which is related to < u 2 >, is 
essentially the phonon propagator, 



(tu) n ) 2 - UJ 2 



, oj n = 2nirT. 



(8) 



With Lo q already defined we have a self consistent equa- 
tion for the fluctuation in the optical mode, 



\2T I 



Which in its asymptotic forms reduces to, 

Td 3 q 

uj 2 — v + gp + vSq 2 + 3Xr<7 

in the high temperature limit and to 

d 3 q 



(9) 
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y/ujQ — v + gp + v5q 2 + 3Xr<7 
in the low temperature limit. Defining the parameters 



(11) 



A («o-«) j 3A'AA 

A = " and Pc = , 

3A q 



(12) 



which give vicinity to the quantum critical point and the 
critical pressure respectively, we can write 

uj 2 {q) = 3AA(1 +p/p c ) - v5q 2 a 2 /2 + 3X R a. (13) 

In absence of coupling to strain, A completely determines 
the quantum phase diagram, with A = giving the 
quantum critical point. Up to this point the result is 
just the renormalization of the factor A to A(l +p/p c ) 
and it becomes an experimentally controllable parame- 
ter. To find out the dependence of T c on pressure, one 
needs to find out temperature dependence of a at cer- 
tain value of p and then T c can be found out from the 
equation u 2 (q,T c ) = 0. Near quantum critical point a 
non-self consistent calculation, with a temperature de- 
pendent momentum cut off (q ma x ~ T) gives a ~ T 2 
which thereby imply T c ~ (1 +p/p c )i. 



III. LANDAU EXPANSION 

Now using this scaling behavior for T c , we switch over 
to a Landau description for the ordered phase, that is, 
we write the variational Gibbs free energy functional as, 



F(P,e:p,T) = ^(T-T c (e))P 



2 1 -P i + -Ke 2 -pe. (14) 



Here a, A > 0, and K is the stiffness constant for strain 
and p is external hydrostatic pressure. P and e are the 
average values of polarization and strain respectively. In 
this scheme we can proceed up to the tricritical point 
along a second order line, i.e. the coefficient of quar- 
tic term is positive and truncation of free energy func- 
tional at quartic term is sufficient. A similar function 
was used recently by Gehring^ to describe various aspects 
of pressure induced quantum phase transition in itiner- 
ant magnets. Starting from equation ([T]), if we consider 
only mean field behavior of both the order parameter and 
strain, we get < e(q) >= j?5{q) and T c ~ (1 +p/p c ). It 
is to be noted that the behavior near QCP can be un- 
derstood only beyond the Landau scheme and a proper 
account of local quantum fluctuations as well as the spa- 
tial variations of order parameter in the free energy func- 
tional is needed. With introduction of one more varia- 
tional parameter (<r), the mean square fluctuations de- 
fined in the equation ([7]) and elimination of it through 
minimization of free energy, will lead to the self consis- 
tent equation in the disordered phase. To a good approx- 
imation we neglect fluctuations in strain (when K 3> 0) 
and considering it's mean field value %8{q), it leads to 
T c ~ (1 +p/p c )?. The reason for such behavior in the 
ordered phase is that in the self consistent calculation 
there a will be replaced by a+ < u > 2 in the expression 
for renormalized frequency ([6|). Here < u > is the mean 
average displacement and the polarization (P) is propor- 
tional to it. Classically < u >^ (T c (p))5 j hence its con- 
tribution retains the form of equation (|13[) same except 
modifying the numerical factor in the right hand side. In 
quantum domain there is possibility that < u >~ T". 
As long as v > i, this assumption is valid. One can also 
question the validity of the linear expansion in temper- 
ature around T c when T c — > 0. In a suitable descrip- 
tion for the phase transition at finite temperature near 
QCP, a higher power of T should appear. Here we fo- 
cus only on QCP at zero temperature. Now we solve the 
coupled equation obtained from the minimization of free 
energy with respective to these parameters which gives 
their physical values. Minimization of free energy with 
respective to polarization gives 



dF(P,e:p,T) 
OP 



= = P [a(T - T c (e)) + AP 2 ] . (15) 



We need one more equation to eliminate e from the above 
equation and this equation is obtained from the mini- 
mization of free energy with respective to strain as fol- 
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lows, 



In terms of pressure this relation becomes 



dF(P,e:p,T) n a dT c (e) „ 



8e 



= ° = -2 



P z + Ke-p (16) 



One needs to solve coupled equation © and © to get 
physically observable quantities. 

Paraelectric Phase: This phase is characterized by 
zero average value of polarization and leads to the so- 
lutions 



P = 0, e=p/K and F {p,T) 



P 

2K 



(17) 



Ferroelectric Phase: In this phase average polarization 
takes non zero value and interplay between strain and 
polarization becomes interesting. Here, 



P 

p_ 
K 



± 



a(T c (e) - T) 



2XK 
5(e) (say) 



A 



(T c (e)-T) 



dT c (e) 
de 



(18) 



and on substituting this expression the free energy in the 
polar phase becomes 



F p ol(p,T) = F (p,T) 



Kg 2 {e ) 



rx (T c (e )-Ty + 2 



F (p,T)-^.(T c (e )-Tf 



a 2 ,dT c (e) 



2XK v de 



(19) 



Where eo is a solution of equation (fl~8"|) . The second term 
in the right hand side of the equation (TT9"|) vanishes at T c 
and also below T c if 
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a 2 ^T c (e), 



2XK K de 



kr = o 



(20) 



Clearly the above equation is the condition for tricritical- 
ity. Thus the condition to achieve first order transition 
is given by 



F pol (p,T) > F (p,T) 



dT c 



2KX V de 



> 1 



(21) 



a 2 fdT c dp 



2KX V dp de 



a 2 K fdT c \ 2 



(22) 



The point in p-T phase diagram satisfying the equality 
marks the transition between second and first order is 
called "tricritical point". In our estimation of T c (e) or 
T c (p) (estimated microscopically and supported experi- 
mentally) is taken as 



T c {p) = {l-p/p f 



(23) 



where £ = i < 1, which tells that at quantum critical 
point the transition is inevitably first order. 



IV. DISCUSSIONS 



The above mentioned analysis is a qualitative answer 
to the recently found phase separation behavior near fer- 
roelectric quantum critical point and its absence when 
one is away from the quantum critical point. Our analy- 
sis shows that the exponent £ which determines the scal- 
ing behavior of T c with pressure near quantum critical 
point, determines whether phase separation is there or 
not. The fact that £ < 1 is clearly a consequence of the 
large quantum fluctuations near QCP. Near QCP system 
becomes very much sensitive to any perturbation and is 
manifested by the divergence of dT c /dp) Pc . These tells 
that in the proposed form of free energy in equation (jT4j) 
some amount of fluctuation effects are included and the 
first orderness and the occurrence of phase separation is 
essentially fluctuation induced. Here we have estimated 
exponent £ in a qualitative manner. Correct estimation 
needs to take higher order correlations into account and 
needs more involved calculations. For the present pur- 
pose only relevant aspect about £ is, whether it is < 1 or 
not. As long as one have £ < 1 (from theory or experi- 
ments) the whole discussion is true. 
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